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Abstract 

We consider the impedance of a structure with rectangular, periodic corrugations on two opposing 
sides of a rectangular beam tube. Using the method of field matching, we find the modes in such 
a structure. We then limit ourselves to the the case of small corrugations, but where the depth 
of corrugation is not small compared to the period. For such a structure we generate analytical 
approximate solutions for the wave number k, group velocity Vg, and loss factor k for the lowest 
(the dominant) mode which, when compared with the results of the complete numerical solution, 
agreed well. We find: if w ~ a, where w is the beam pipe width and a is the beam pipe half- 
height, then one mode dominates the impedance, with k l/\/^w5 {5 is the depth of corrugation), 
{1 — Vg/c) ~ 5, and n ~ l/^aw), which (when replacing w by a) is the same scaling as was found for 
small corrugations in a round beam pipe. Our results disagree in an important way with a recent 
paper of Mostacci et al. [A. Mostacci et al, Phys. Rev. ST-AB, 5, 044401 (2002)], where, for the 
rectangular structure, the authors obtained a synchronous mode with the same frequency k, but 
with K ^ 5. Finally, we find that if w is large compared to a then many nearby modes contribute 
to the impedance, resulting in a wakefield that Landau damps. 
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I. INTRODUCTION 



In accelerators with very short bunches, such as is envisioned in the undulator region 
of the Linac Coherent Light Source (LCLS) the wakefield due to the roughness of the 
beam-tube walls can have important implications on the required smoothness and minimum 
radius allowed for the beam tube. One model that has been used to study roughness is a 
cylindrically-symmetric structure with small, rectangular, periodic corrugations. For such a 
structure, if the depth-to-period ratio of the corrugations is not small compared to 1, it has 
been found that the impedance is dominated by a single strong mode with wave number 
k ~ l/^/a6, with a the structure radius and 6 the depth of corrugation, and loss factor 
K = (in Gaussian units) [|, |]. 

In a recent report Mostacci et al. [Q, studied the impedance of a structure with small, 
rectangular, periodic corrugations on opposing sides of a rectangular beam tube using a 
perturbation approach. For a beam tube with width w comparable to height 2a the authors 
find a mode with a similar frequency dependence as in the round case, but with a loss factor 
that is proportional to the depth of corrugation 5. If this model is meant to represent surface 
roughness with e.g. 5 ~ 1 /im and a ~ 1 cm, then their result implies a factor ~ 10~^ smaller 
interaction strength than was obtained in the earlier cylindrically symmetric calculations. 
Such a result seems unlikely — we would not expect a huge difference in loss factor when 
changing from round to rectangular geometry. It is the goal of this paper to resolve this 
discrepancy and to show that a correct calculation for the rectangular cross section indeed 
gives a result that differs only by a numerical factor from the round case. 

Another motivation for this work is to understand the impedance of two corrugated plates, 
the limit of our geometry when w becomes large. And although, when w is not large, the 
geometry is somewhat artificial, it may still be a useful model for some vacuum chamber 
objects of accelerators, e.g. for the screens in the LHC vacuum chamber [Q. And thirdly, 
we note that fabricating a structure with artificially large corrugations, for the purpose of 
experimentally studying roughness impedance, may be much easier for the rectangular than 
the round beam pipe. 

In this report we calculate the impedance of the rectangular structure of Mostacci et 
al. — but not limiting ourselves to small corrugations — using the method of field matching. 
The solution is written as an infinite homogeneous matrix equation that we truncate to 



2 



solve numerically. Note that our approach is very similar to that used for the analogous 
cylindrically symmetric problem in the computer program TRANS VRS [Q. Note also that 
recently, Xiao et al. used a similar method to solve the impedance of the rectangular 
structure, but with the corrugated surfaces replaced by dielectric slabs 0. Next, using a 
perturbation approach applied to the field matching equations we find the analytical solution 
for the limit of small corrugations. Finally, we compare the analytical to the numerical 
results. 



II. FIELD MATCHING 



We consider a periodic, rectangular structure with perfectly conducting walls, two periods 
of which are sketched in Fig. [1|. In the horizontal [x) direction the structure does not 
vary, except for walls at x = ±w/2. One period of the structure extends longitudinally to 
z = ±p/2. This cell can be divided into two regions: Region I, the "tube region", extends 
to y = ±a; Region II, the "cavity region", for z = ±g/2, extends beyond y = ±a to 
y = ±(a + (5). An exciting point beam moves at the speed of light c from minus to plus 
infinity along the z axis. We are interested in the steady-state fields excited by the beam, 
and assume that initial transients have all died down. Note that we will work in Gaussian 
units throughout. 
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FIG. 1: A longitudinal cut of the structure geometry considered here, showing two periods in the 
z-y plane (left), and a transverse cut showing the cross-section of the structure (right). 



We assume that the fields of a mode excited by the beam have a time dependence e-^'^'^*, 
where k is the mode wave number and t is time. For either region the fields can be obtained 
from two Hertz vectors, 11^ and He, which generate, respectively, TM and TE components 
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of the fields: 



£ = VxVxn^- jkV X He , (1) 

u = vxVxne + jkv X . 

Since there is no variation in the x direction we choose it as the direction of the Hertz 
vectors. To satisfy the boundary conditions aX, x — ±u'/2 the fields vary as cosines and sines 
of kxX where 

= - , (2) 

w 

with m an odd integer (see below). The general solution involves a summation, over all m, 
of such modes. 

Consider modes with horizontal mode number m. In the tube region, the most general 
form of the (x component of the) Hertz vectors, consistent with the (perfectly conducting) 
walls aA, X = ±w/2, and the Floquet condition in z is: 

oo 

nL = E Ksinh(fc;„|/) + i?„cosh(fc;„|/)]sin(/c,a;)e-^-^"^ , (3) 

n=— oo 

oo 

nL = E [Cn sinh(A;i„|/) +Z}„cosh(/c;„l/)] cos(/c^a;)e-^'^"^ , 

n=— oo 

with 

/3n = /?o + ^ , kln = ^(^l-k^ + kl . (4) 

Since the structure is symmetric in y about y — 0, the field components will be either even 
or odd in y, and the modes will split into two categories. In the first type — — and 
the resulting modes have Sz on axis, in the second type Bn — Cn — and the resulting 
modes have = on axis. In either case we are left with only 2 sets of unknown constants 
in Region I. Since an on-axis beam can only excite modes of the first type, it is this type in 
which we are interested. 

In the cavity region, the most general form of the Hertz potentials, consistent with per- 
fectly conducting boundary conditions at z — ±gi/2 and y — ±(a + 5) is: 

oo 

^mx = E^^si^K(« + <^-y)]sin(^^^)sin[Q;s(^ + ^/2)] , (5) 

s=l 
oo 

^ii = J^Fs cos[kll{a + 5-y)] cos{k^x) cos[as{z + g/2)] , 



s=0 
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with 

a, = y , kll = ^k'^ -al-kl . (6) 

Note that in both regions £y, E^, and Tij. depend on x as cos(A;^x) and, therefore, the 
boundary conditions on the walls at x = are automatically satisfied. 

We need to match the tangential electric and magnetic fields in the matching planes, at 
y = ±a: 

^ (7) 

[ : g/2<\z\<p/2 

= nil, ■■ 1^1 < ^?/2 • (8) 

Using the orthogonality of e~^^"^ over [— p/2,j9/2] in Region I, and sin[as(z + g/2)\ and 
cos[a;s(z + 9/2)] over [—g/2,g/2\ in Region II, we obtain a matrix system that we truncate 
to dimension 2(2A/' + 1) x 2(2A/' +1), where M is the largest value of n that is kept. To 
obtain modes excited by the beam we need to set (3n = k for one value of n. The frequencies 
at which the determinant of the resulting matrix vanishes are the excited frequencies of the 
structure. 

The relation of the coefficients at the excited frequencies gives the eigenf unctions of the 
modes, from which we can then obtain the (i?/Q)'s and the loss factors. The loss factor, 
the amount of energy lost to a mode per unit charge per unit length of structure, is given 
by 



|2 



K 



zs 



(9) 



Au{l-vjc) ' 

with £zs the synchronous component of the longitudinal field on axis, u = 
(87rp)~^/ \S\'^dx dy dz, the (per unit length) stored energy in the mode [the integral is over 
the volume of one period of structure], and Vg the group velocity in the mode. Note that 
the factor 1/(1 — Vg/c) is often neglected in loss factor calculations (it appears to have been 
neglected in Mostacci et al). This factor in the loss factor, which — as we will see — is very 
important in structures with small corrugations, is discussed in Refs. 0, ||, ^J; we give a new 
derivation of it in Appendix A. Finally the longitudinal wakefield is given as 

I^(s) =2$(s)^fi:„cos(M , (10) 

n 

with $(s) = for s < 0, 1 for s > 0, and the sum is over all excited modes. 
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In Appendix B we present more details of the calculation of the modes of the corrugated 
structure using field matching. We have written a Mathematica program that numerically 
solves these equations for arbitrary corrugation size. The results of this program will be 
used to compare with small corrugation approximations presented in the following section. 

III. SMALL CORRUGATIONS 

Let us consider the case where the corrugations are small, but with 5~5f~p<Ca~iy. In 
the analogous cylindrically symmetric structure it was found that: (i) there is one dominant 
mode (its loss factor is much larger than those of the other modes), (ii) this mode has 
a low phase advance per cell, and (iii) the frequency of the mode k ~ l/^/a5 [|], p!0| . 
For our rectangular structure we look for a mode with the same properties. As was the 
case for the cylindrically symmetric problem we also assume that the fields in the cavity 
region are approximately independent of z, and that one term in the expansion of the 11 
vectors, the term with n = and s = 0, suffices to give a consistent solution to the field 
matching equations [Q. Note that, it is true that to match the tangential fields well on the 
matching plane may require many space harmonics (though even then, near the corners, 
Gibbs phenomena and the edge condition will result in poor convergence); nevertheless, as 
with the analogous cylindrically symmetric problem, the global mode parameters in which 
we are most interested — frequency k, group velocity Vg, and loss factor k — can be obtained 
to good approximation when keeping only the one (the n = 0, s = 0) term. 

Setting a = implies that 11^^ = 0, and that there are only 3 non-zero field components 
in the cavity region: S^^ , Sy^ , and Ti^J ■ For small corrugations the excited modes become 
approximately TM modes. To allow matching at the interface of Regions I and II we end 
up with 

nL ^ , (11) 

^ Co smh{kl,y) cos{k,x) e"^'^"^ , 

and 

^ , (12) 
Ilg^ ^ jFo cos[fc^o(a + 6 -y)] cos^k^x) . 
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Let us sketch how we match the fields: We equate 8z and Tix for the two regions at 
y = ±a; we multiply the first equation by e^^°^ and integrate over one period in z, and then 
we integrate the second equation over the gap in z. When we divide the resulting equations 
one by the other, the constants Cq, Fq, drop out, and we are left with an approximation to 
the dispersion relation, one valid in the vicinity of the synchronous point (the subscript 
for P is understood): 

y^^^^fc^T^ coth + = ^^^^^^^yp^tan (^^k^-kl6^ . (13) 

To properly keep track of the relative size of the terms in further calculations, we assign 
to each parameter an order using the small parameter e: let a, w, be of order 1; 6, g, p, of 
order e^; and k, (3, of order 1/e. To find the synchronous frequency we let (3 = k in Eq. (jlSj), 
expand the equation to lowest order in e, and then set e = 1. The result is 

kl, = ^ coth (k^a) , (14) 

(the subscript m is included here to remind us of the m dependence). Note that, if a ~ 
(and p ^ g) then k ~ l/\fw5^ which is of the same order as the result that was found for the 
cylindrically symmetric problem. Note also that, for the limit g = p, the dispersion relation 
and the synchronous frequency here agree with those given in Mostacci et al. 

For the group velocity we take the partial derivative of Eq. (|13D with respect to /5, and 
rearrange terms to obtain 1 — dk/d(3 = 1 — Vg/c. After expanding in e, keeping the lowest 
order term, and finally setting e = 1 we obtain 



^ _ {vg)m \ _ 2(5 k^g 



c I p 



sinh^(/Ca:a) 



(15) 



?,mh.{kxa) cosh(fc^a) — k^a 
Note that, as in the cylindrically symmetric problem, (1 — Vg/c) ~ 5. The loss factor of our 
structure 

27r 

i^m = — F (k^a) , (16) 
wa 

with 

Fix) = . u ^ u ^ ■ (1^) 

smh(xj cosh(xj 

The function F{x) and an approximation for large x shown in Fig. |^. Note that for 
K in the MKS units of [V/pC/m], one multiplies Eq. (|16D by the quantity Zoc/(47r), with 
Zq = 377 Q. Note also that our result is independent of 6, unlike the result of Mostacci et 
al. 
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FIG. 2: The function F{x) (solid) and the approximation 4xe~^-^, vahd for x ~ 1 (dashes). 

The total longitudinal wakefield is given by Eq. (|r^). Note that, if w ^ a (% > 1) then 
one mode dominates the wake, just like in the round case. (For example, if x = 1, then the 
amplitude of the first, m = 1 term is 20 times larger than that of the next, m = 3 term 
in the wake sum.) If, however, w ^ a, then more than one mode will contribute to the 
impedance of the structure; in the limit of w oo (two corrugated plates) there will be a 
continuum of modes contributing to the impedance. The impedance is given by the Fourier 
transform of the wake. Its real part is 

TZeZ = TT ^ Km — kjnc) + S{uj + kmc)] . (18) 

m 

Consider now the limit of two corrugated plates {w —* oo). The mode spectrum becomes 

continuous and the sum in Eq. (|TBp can be replaced by an integral 

Tc r°° 
neZ = - dxFix) 
Jo 

The integral can be solved numerically, with the use of the relation / dx g{x)6[f{x)] = 
[g{x)/\f'{x)\]x=xo where /(xq) = 0. The result is shown in Fig. |a; note that the axes are 
normalized to kr = \Jp/ (aSg) and = n/^a^krc). We see a continuous spectrum of modes 
beginning at wave number kj., with average l.lAkr and rms 0.18/cr- The corresponding 
wakefield becomes a damped oscillation (see Fig. Qd). We see an effective Q ~ 10. Note 
that W{0^) = 7r^/(4a^) [to be discussed more in a later section]. 

Finally, we should point out that it has been observed for the case of the cylindrically 
symmetric problem that, if the small corrugations are replaced by a thin dielectric layer of 
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kr = \/p/JaSg) and Zr = i: / {a?krc). 



thickness 5, and if the correspondence is made that the dielectric constant e = p/{p — g), 
then the resuhs for the two problems are the same P]. Recently the modes in a rectangular 
structure of Fig. |l|, but with the corrugated surfaces replaced by dielectric slabs, have been 
obtained by Xiao et ai, also using a field matching approach 0. If we take their results, 
letting the thickness of the dielectric layers (5) be small, we obtain our results for k, Vg, and 
K when we make the correspondence e = p/{p — g)- 

A. Comparison with Numerical Results 

To test the validity of the analytical approximations in the case of small corrugations, 
we compare with numerical results obtained by the Mathematica field matching program 
(the method of solution is described in Appendix B). Consider as an example a square beam 
tube {w/a = 2) with p/a = 0.05, g/a = 0.025, and 6/a = 0.025, and let us consider the 
lowest (m = 1) mode. In the field matching program we take S = A and N" = 4, i.e. 5 
space harmonics are kept in the cavity region and 9 in the tube region. (We find that, for 
the example geometry, keeping more terms has no significant effect on the results.) 

We begin by comparing the dispersion curve (see Fig. ^). Shown are the field matching 
result (the solid curve) and the approximation, Eq. ([T3|) (the dashes). We see that the two 
agree well except far from the synchronous phase. The cross plotting symbol locates the 
synchronous point, with kp = 0.2007?, a result which is 7.5% larger than the analytical value 
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of Eq. (^4|). It is interesting to note that this dispersion curve is almost identical to the one 
obtained (also by field matching) for the same geometry but in a round beam pipe |^ . As 
for the loss factor, we find that it is a factor 0.84 as large as the analytical approximation, 
Eq. 0. 
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/Sp/TT 

FIG. 4: A dispersion curve example: shown are the numerical result (solid), the synchronous 
point (the cross plotting symbol), and the approximation, Eq. ( p!3|) (dashes). Also shown is the 
speed of light line (dots) . 

These results confirm the validity of the analytical approximations for the structure with 
small corrugations, provided that the depth of corrugation 6 is not small compared to the 
corrugation period p. However, in Ref. it was shown that for the analogous round 
structure the corresponding analytical formulas break down when 6 becomes small compared 
to p: as 6 decreases the frequency first increases than decreases as compared to the analytical 
result; meanwhile the loss factor continually decreases. When 6 is small compared to p the 
impedance is no longer well characterized by a single resonance, and is best described by 
a different model [^. As expected, we find the same kind of behavior in our rectangular 
structure. If, for example, we reduce 6 in our example problem by a factor of 2, we find that 
the frequency becomes 18% larger, and the loss factor 30% smaller, than the values given 
by the analytical formulas. 
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B. Discussion 



Our result for the loss factor, Eq. ([T6|), is independent of the depth of corrugation 5, as 



was found previously for the analogous cylindrically symmetric problem IC]. This result, 
however, is in disagreement with the result of Mostacci et ai, where the loss factor was 
found to be directly proportional to 6. This discrepancy is important to resolve. 

There is a general relation that holds for the wake directly behind the driving particle 

2 

W{0+) = - 7^eZ(cJ)rfcu = 2 Vk„ , (20) 
n Jo 

a relation that does not depend on the specific boundary conditions at the wall. To discuss 
it, consider first the analogous cylindrically symmetric problem. It was earlier found that, 
as long as the corrugations are small and the depth 6 ^ p, the contribution of one mode 
dominates the wake sum. In this case, it was found that, as here, ly (0^) (or n) is independent 
of (5 0. If the corrugations are replaced by a thin dielectric layer, 14^(0"'') does not depend 
on the dielectric properties (neither 6 nor e) [T^ . In the same way, if the corrugations are 
replaced by a lossy metal, 1^(0+) will not depend on the conductivity |13|. And in all three 
cases the answer is the same: iy(0"^) = [In fact, this relation is also valid for the 



(steady-state) wake of a periodic accelerator structure, with a the iris radius ||T4|, |15 



We expect the same type of behavior to hold in a corrugated, rectangular structure, 
i.e. that 1^(0"'") depends only on the cross-section geometry of the beam pipe. In Fig. |^ 
we plot, for our rectangular structure, a^W^(0+)/2 = a^J^m^m, as function of na/w (the 
solid curve). Also shown is the contribution of only the first (m = 1) term (dashes), and 
the approximation 8(7ra/w)^ exp(— 27ra/u') (dots). Note that, for ira/w small, many modes 
contribute to the sum; for na/w ^ 1, one mode dominates. As with the cylindrically- 
symmetric case, 1^(0"'") must still be correct if we replace the corrugated surfaces by thin 
dielectric slabs, or by lossy metal plates. We know of no published result for 1^(0"'") in 
our rectangular geometry to compare with; nevertheless, Henke and Napoly found 1^(0"*") 
between two resistive parallel plates [|16|, which becomes the limit of our geometry as w — > oo. 
Their result, a^W {0^)/2 = 7r^/8, agrees with our calculation for na/w — * 0, and confirms 
our result. 

For a given bunch shape and fixed 6/p, as the depth of corrugation 6 decreases, we expect 
the induced voltage (the convolution of the bunch shape with the wake) to also decrease. 
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FIG. 5: The sum of the loss factors a'^J2m'^m [= a^VF(0^)/2] as function of ira/w (sohd). 
Also shown are the contribution of the first mode, o^ki (dashes), and the approximation 
8(7^a/^^;)^ exp(— 27ra/u;) (dots). 

If the loss factor does not depend on 6 how does this happen? The answer is that as 6 
decreases, the mode frequency k increases, and the wake, when convolved with the bunch 
shape, will yield an induced voltage that will decrease (at least as fast as l/k). Concerning 
this question, the wake of this structure behaves similarly to the resistive wall wake for very 
short bunches as the conductivity increases: 1^(0"'") also does not change but the wake first 
zero crossing moves closer to s = 0. 

IV. CONCLUSION 

We studied the impedance of a structure with rectangular, periodic corrugations on two 
opposing sides of a rectangular beam tube using the method of field matching. We described 
a formalism that, for arbitrary corrugation size, can find the resonant frequencies k, group 
velocities Vg, and loss factors k. In addition, for the case of small corrugations, but where 
the depth of corrugation is not small compared to the period, we generated analytical per- 
turbation solutions for k, Vg, and k, for the dominant mode. We then compared, for such a 
structure, the results of the computer program and the analytical formulas, and found good 
agreement. 

In general, we found that, for the structure of interest, the results are very similar to what 
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was found earlier for a structure consisting of small corrugations on a round beam pipe: if 
w ^ a, where w is the beam pipe width and a is the beam pipe half- height, then one mode 
dominates the impedance, with k ~ l/^/a6 {6 is the depth of corrugation), (1 — Vg/c) ~ 5, 
and K, ~ If, however, w is large compared to a we find that many nearby modes 

contribute to the impedance, resulting in a wakefield that Landau damps. 



APPENDIX A: EXCITATION OF A SYNCHRONOUS MODE BY A MOVING 
RELATIVISTIC POINT CHARGE 

Consider first a cavity of frequency u with the electric field of an eigenmode S{r)e^^^. 
The energy in the eigenmode is denoted by U . If a point charge q passes through the cavity, 
it excites this mode to the amplitude As (where As is a complex number), so that after the 
passage through the cavity the electric field of the mode will be AsS{r)e^'^^ , and the energy 
lost by the charge is equal to lA^pf/. In quantum language, this is spontaneous radiation 
of the charge into the mode under consideration which is indicated by the subscript s. It is 
clear that As is proportional to the charge of the particle q. 

To calculate the amplitude As, let us consider a situation when, before the charge enters 
the cavity, the latter already has this mode excited by an external agent (RF source) to the 
amplitude Aq. Due to linearity of Maxwell's equation, after the passage of the charge, the 
field in the cavity will be equal to the sum of the initial mode Aq and the spontaneously 
radiated mode As-, with the energy given by + Aopf/. The change of the energy in 
the cavity is 

Aiy = \As + Aopf/ - |v4opf/ = {AsAq + c.c.)U + \As\'^U , (Al) 

where c.c. denoted a complex conjugate. Let us consider the limit of small charges, g — >^ 0, 
then we can neglect the last term on the right hand side of Eq. (|A1|) , which scales as g^, 
and keep only the first term that is linear in g, 

= {AsAq + c.c.)U . (A2) 

Discarding the term oc q^ means that we neglect the beam loading effect. 

We can now balance the energy change /SW with the work done by the external field Aq 
during the passage of the charge. This work is equalt to the integral of the electric field 
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£z{z) along the particle's orbit 

AW = -qUeAo J dz£,{z)e''''/'' = J dzf,(z)e^""/" + c.c. (A3) 

Comparing Eq. (^) with Eq. ( |A3| ) we conclude that 

^s = -^j dz£z{z)e^-^/\ (A4) 

Hence we found the amplitude of spontaneous radiation of the particle in terms of the 
integral along the particle's orbit of the electric field. 
The energy lost by the particle (loss factor) is 

m^U^qi, (A5, 

where the voltage V = J dz£z{z)e^'^^/'" . 

Let us now apply the same approach as above to the excitation of a mode that propagates 
with the speed of light in a waveguide. To deal with a mode of finite energy we consider 
a wave packet, and assume that the packet has a length L, as shown in Fig. ^ below. It 

L 



FIG. 6: The shape of the wave packet of the synchronous mode. The packet has a long plateau of 
length L and short edges. 

propagates in the pipe with the group velocity Vg. The energy in the mode U can be related 
to the energy flow P (integrated over the cross section averaged over time the Pointing 
vector) if we note that U/L is the energy per unit length, and hence {U/L)vg is the energy 
flow equal to P, hence 

U = —. (A6) 

Now, the particle is synchronous with the wave and stays all the time in the same phase, 
so it sees the same longitudinal electric fleld Sz which we denote by Szs- The integral from 
Eq. ( |A4| ) can be written as 

JdzSz{z)e''^'/'' -^cTSzs, (A7) 
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where T is the interaction time between the wave and the particle. This is actually the time 
when the particle stays in the wave, and taking into account that the wave is moving with 
velocity Vg and the particle is moving with c 

T = — (A8) 

c-Vg 

Hence, for the amplitude of the radiated wave we find 

As = -:^S.s^, (A9) 

2U C — Vg 

and the energy W radiated by the particle 

W = AlU = ^Sl-^^. (AlO) 

To find the energy radiated per unit length of the path, we divide W by the length of 
the interaction path Lc/{c — Vg), which gives 

= |lfi-^, (All) 

where the energy per unit length of the path u = U/L. Finally, since the loss factor 
K = q~'^dW/dz, we arrive at Eq. @. 



APPENDIX B: FIELD MATCHING, THE GENERAL SOLUTION 

In Section II we presented Hertz vectors and wave numbers for Regions I and II, and 
also the four equations that need to be matched at the interface y = ±a. We continue with 
the notation introduced there: We multiply the matching equations for Sz and by e^^"'^ 
and integrate over [— p/2,p/2]; and we multiply the matching equations for Ti^ and by 
sm[as'{z + g/2)] and cos[as>{z + g/2)] and integrate over [—g/2, g/2]. We obtain the infinite 
set of equations: 

(C;fcfc>-i?>7r/5„)cosh(fcJ„6) = ^J2Nnsi-Fikkliw + E'^m7ra,)smikliS) 

P s 
P s 

{E'^kkliw + Flm-nas) cos{kli5) = 2 ^ M,n{B'^kkl^w - C„m7r/3„) sinh(A;^„6) 

n 

(l + 5,o)i^;cos(A;^,^5) = -2 5:c;Ar,„sinh(A;J„6) . (Bl) 

n 
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Here 



E F 









2 




hi 







sin(/5„5(/2) 
cos(/3„5f/2) 



: s even 
: s odd 

s even 
s odd 



(B2) 



(B3) 



and 5,,/ the Kronecker delta. 



This system of equations can be written as a homogenous matrix equation: 



G{H^-P) -GH 
-GH G 




P{Q^ + R^)/S -PQ 
-PQ/S P 





with superscript T indicating the transpose of a matrix. The diagonal elements of di- 
agonal matrices are: Gn — coth{ky,^fi) / (kky^w) , if„ = irnrPn, In — kky^w; Pg — 
2gtan{kyl6)/{pkkylw), Qs = mnas, Rs = kkylw, Ss = {I + Sso)- Note that the sys- 
tem matrix is real. The expansion coefficients are: B'^ = — smh.{ky^b)C'^ and C'^ — 
sinh{klJ)){kkl^wB'^ - m7r/5„C;). 

To solve the matrix equation we truncate to dimension 2(2jV-|- 1) x 2{2J\f + 1), where N 
is the largest value of n that is kept. Therefore, subscript n, representing space harmonic 
number in the tube region, runs from —H to H] subscript s, representing space harmonic 
number in the cavity region, runs from to 5, the largest value kept. Note that the values 
A/", iS, should be chosen so that (2A/'+ ~ (5 + I)/'?- The system matrix f/ is a function 
of I3q and of k. To find synchronous modes, we need to first set, for one space harmonic 
n', /?„/ = k and then numerically search for the value of k for which the determinant of U 
becomes zero. The value n' should be taken to be the nearest integer to kp/ {27r). To find 
values of the dispersion curve, we, for various values of /3„/ [where again n' is the nearest 
integer to kp/(27i)], numerically search for the value of k for which the determinant of U 
becomes zero. 

Once we have found the frequency we can find the eigenfunctions, from which we obtain 



ISzsl"^ on axis. 



I ^?:.<i I — k 



(B5) 



(where ns represents the synchronous space harmonic) and the energy per unit length u. 
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For example, the stored energy in Region I is given by 
1 



u 



[-2k'kla + (k'kl^ + k'(3t) sinh(2fcj„a)/fcij - AB'^C'^kkA smh(2fc;„a)) (B6) 

with k'^_ — k^ — kl, with a corresponding equation giving the energy stored in Region II. 
Note that for small corrugations, u^^ <S . The quantity 1/(1 — Vg/c) is obtained by 
first calculating the dispersion curve, and then finding the slope at the synchronous point 
numerically. Knowing l^^^p, u, and 1/(1 — Vg/c) we can finally obtain the loss factor k. 
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